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DISCRETE  ALGORITHM  FOR  CONTROLLING  THE  FINAL  DESCENT  VELOCITY  OF 
SPACE  VEHICLES  IN  THE  ATMOSPHERE  OF  MARS 


I.  K.  Bazhinov,  N.  M.  Ivanov,  A.  I.  Martynov 


I 


An  analysis  is  made  of  a discrete  algorithm  for 
controlling  the  final  descent  velocity  of  space 
vehicles  (KA)  [HA]  in  the  atmosphere  of  Mars. 

This  can  be  realized  by  simple  automatic  devices. 
Control  of  the  vector  of  lift  is  carried  out  by 
means  of  changing  the  roll  attitude  (i.e.,  the 
effective  component  of  lift).  The  control  algo- 
rithm makes  use  of  the  lines  of  change-over  which 
are  stored  in  the  on-board  computer.  Numerical 
results  are  presented  from  the  calculation  of  the 
effectiveness  of  the  control  algorithm.  It  is 
shown  that  the  control  algorithm  can  be  used  in 
the  construction  of  systems  for  control  of  the 
final  descent  velocity  for  a wide  class  of  landers 
(SA)  [CA]  and  for  different  entry  velocities. 


In  work  [1]  data  were  presented  concerning  the  optimal  control 
of  the  final  descent  velocity  of  a space  vehicle  in  the  atmosphere 
of  Mars.  In  particular  it  was  shown  that  for  a certain  class  of 
landers  the  optimal  control  program  (or  close  to  it)  is  a program 
of  single  change-over  of  "effective”  quality  cosy).  This 

program  is  realized  quite  readily  bv  simple  means  [2].  However, 
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the  program  of  control  with  a single  change-over  of  has  limited 
application,  since  it  can  be  used  only  for  a comparatively  narrow 
class  of  landers  with  low  values  of  reduced  load  on  the  front  sur- 
face Px^250  kgf/m2  and  magnitude  of  lift-drag  ratio  Ky^0.3. 

Figure  1 shows  the  dependence  of  final  velocity  on  load  on  the 
frontal  surface  in  the  case  of  different  programs  for  control  of 
roll  attitude.  Along  with  this,  the  investigations  showed  that 
for  a wide  class  of  landers  (px^250  kgf/m2  and  Kf^.0.3)  the  mini- 
mal value  of  final  velocity  is  realized  when  using  a simple  program 
of  double  change-over  of  "effective”  quality.  Below  we  will  con- 
sider a discrete  algorithm  of  control  which  incorporates  a twofold 
transfer  of  quality  K . . 

3<P 

We  will  accept  that  the  means  of  navigation  ensure  an  accuracy 

of  guidance  which  is  sufficient  for  entry  into  an  operational 

approach  corridor  [1],  Just  as  previously  [l]-[3],  in  all  cases 

we  will  assume  that  the  minimal  permissible  altitude  of  flight  of 

the  lander  above  the  surface  of  Mars  H . is  no  less  than  the 

min  flon 

altitude  of  the  beginning  of  operation  of  the  soft  landing  system: 

Hm\tt  im  ^ Ha.  ^ ^ 

When  using  a program  of  twofold  transfer  of  quality  a lander 
enters  the  atmosphere  with  Kg^+K^.  (y=0),  after  a certain  time 
follows  the  transfer  Kg^-K^  and  then  again  to  K^^+K^.  The 
moments  of  transfer  are  selected  from  the  condition  of  obtaining 
V H min  with  the  simultaneous  fulfillment  of  limitation  (1). 


2. 


r 

VM{*/cex] 


SU  pxpte/M 

(P 


Figure  1 


Key:  (1)  V [m/s];  (2)  V = 6 km/s;  H =6  km;  K*=0.3 ; 

K BX  K w 

(3)  optimal  control;  (4)  single  change-over  y ; (5)  twofold 

p 

change-over  y ; (6)  Px[kgf/m  ]. 

For  realizing  the  program  of  twofold  change-over  it  is  neces- 
sary in  each  specific  case  to  be  able  to  determine  the  two  points 
of  change-over  by  using  on-board  devices.  This  can  be  realized  by 
using  two  lines  of  change-over  which  are  stored  on  board  the  lander, 
The  method  of  selection  of  the  lines  of  change-over  is  analogous 
to  that  described  previously  for  the  case  of  a single  change-over 
of  quality  [2].  As  an  example  we  will  consider  the  algorithm  for 
twofold  transfer  for  controlling  the  final  descent  velocity  V 

H 

o 

of  a lander  with  Px=350  kgf/m  and  Kf  =0.3.  The  lines  of  change- 
over for  this  case  are  shown  in  Figure  2 in  the  form  of  dependences 
of  acceleration  nTvt(V„_)  and  apparent  velocity  t_(V  _)  [2].  It  is 
evident  that  the  lines  of  change-over  t (Vs  ) have  a simpler  form 
than  nxn(vsn)>  and  accordingly  they  can  be  approximated  more  easily 


by  polynomials  of  the  1st  and  2nd  degree.  Considering  this  fact, 
and  also  the  fact  that  the  method  of  assigning  the  lines  of  change- 
over has  virtually  no  influence  on  the  magnitude  of  V [2],  subse- 

H 

quently  we  will  assume  that  the  lines  of  change-over  on  board  the 
lander  are  assigned  in  the  form  tn(Vsn). 

For  this  particular  case  the  lines  of  change-over  shown  in 
Figure  2 can  be  approximated  in  the  following  manner.  The  line  of 
change-over,  determining  the  moment  of  the  first  transfer  - by  one 
straight  line: 

bn  1 = al  Vsn  + bl» 
where  a-^  = -11,111  s2/m,  b1  = 101,666  s. 

The  line  of  change-over,  determining  the  moment  of  the  second 
transfer  - by  three  straight  lines: 

*n  21  * a21  vs„  + »21  <V,<V8l). 
where  a2l  * -113,461  s2/m,  b2l  = 481,942  s,  Vsl  = 3.45  km/s; 

n 22  = a22  Vsn  + b22  ^sl^s^s2)t 
where  a22  = -925  s2/m,  b22  = 3191,25  s,  Vg2  = 3.47  km/s; 

tn  23  = a23  Vsn  + b23  ^Vs>Vs2^» 

2 

where  a2^  ■ 100  s /m,  b2^  = - 275  s. 


The  length  of  the  time  interval  At  between  the  two  moments  of 
change-over  k^  for  each  type  of  lander  depends  significantly  on 
the  altitude  of  the  conditional  pericenter  of  entry  trajectory  . 
Thus,  for  example,  for  the  lander  being  considered  in  the  case  of 
movement  near  the  upper  boundary  of  the  corridor  At  has  the  great- 
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est  value  and  comprises  approximately  70  s.  A reduction  of  leads 
correspondingly  to  a lessening  of  the  time  Interval  between  the 
two  transfers  right  up  to  the  complete  disappearance  of  it  dur- 
ing movement  of  the  lander  on  the  lower  boundary  of  the  approach 
corridor  H^.  In  the  last  case  the  craft  is  moving  on  the  entire 
descent  trajectory  with  an  unchanged  value  Kg^  = + . 


Figure  2. 

Key:  (1)  px=350  kgf/m2;  Kf=0.3i  VBX=5.8  km/s;  (2)  tn  [s] ; 

(3)  Vsn  [km/s]. 

The  table  contains  the  results  of  an  appraisal  of  the  effect- 
iveness of  functioning  of  the  algorithm.  Also  presented  there  are 
the  maximum  magnitudes  of  the  main  perturbations. 
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Table 


Px=350  kgf/m2;  Ymax=";  Kf  =0.3;  Hh  = 6 km;  = -70  km; 

V =5.8  km/s. 

B X 


Key:  (1)  Type  of  perturbation;  (2)  Range  of  perturbation; 

(3)  VK,  [m/s]. 

It  is  evident  that  the  use  of  the  proposed  algorithm  makes  it 
possible  to  work  out  all  the  perturbations  quite  well.  Mathemati- 
cal expectation  MCV,,)  and  triple  mean  square  deviation  of  the 
final  velocity  aVh=3 o,  determined  with  the  possible  combined  action 
of  all  the  perturbations  according  to  the  method  of  B.  G.  Dostupov 
[1],  [4] , comprise  correspondingly  M(VU)=638  m/s  and  AV  =45  m/s. 

The  possible  error  in  the  operation  of  the  algorithm  for  con- 
trol of  final  velocity  does  not  exceed  60  m/s  In  comparison  with 
"ideal"  control  (instantaneous  transfer  of  the  lander  for  bank, 
absence  of  Instrument  errors,  etc.),  and  120  m/s  in  comparison  with 
optimal  control  (see  Figure  1). 

It  is  necessary  to  note  that  the  minimum  flight  altitude  for 
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the  lander  on  a perturbed  trajectory  can  take  values  which  differ 
from  the  assigned  Hmin  n>  Thus  f’or  the  example  being  considered 
the  magnitude  of  mathematical  expectation  of  the  minimum  flight 
altitude  under  the  possible  total  action  of  all  perturbations  com- 
prises M(Hmin)=6-2  km,  and  the  magnitude  of  the  triple  mean  square 
deviation  3a(Hmin)=0- 88  km. 

Px=350  kgf/m*;  =0.3;  Vbx=5.8  km/s . 


Figure  3. 

Key:  (1)  VK[m/s]j  (2)  [rad]. 

Now  let  us  consider  how  the  change  in  the  maximum  magnitude 
of  the  roll  attitude  Ymax  and  the  minimal  permissible  altitude  of 
flight  above  the  surface  of  Mars  influence  the  accuracy  of  per- 
formance of  the  control  algorithm..  As  is  evident  from  Figure  3, 
a reduction  of  Ymax  (in  the  case  of  the  corresponding  change  of  the 
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* 


lines  of  change-over,  see  Figure  2)  leads  to  an  increase  of  M(V  ), 
and  in  this  case  the  magnitude  of  3a  remains  virtually  unchanged. 
Thus,  for  example,  with  a change  of  y_„„  from  tv  to  tt/2  the  magni- 
tude  of  M(V  ) increases  from  638  to  697  m/s,  and  AV  =3a  from  45 

H K 

to  50  m/s.  An  increase  in  the  minimal  permissible  flight  altitude 
leads  to  an  increase  of  M(VH)  and  AV^ . Thus  with  a change  of 
H . „ from  3 to  9 km  (lines  of  change-over  for  these  cases  are 

shown  in  Figure  2)  the  magnitude  of  M(V  ) increases  from  518  to 

K 

757  m/s,  and  AV  from  23  to  58  m/s.  Thus  the  materials  presented 
testify  to  the  sufficiently  high  effectiveness  of  the  algorithm 

p 

with  two  lines  of  change-over  for  landers  with  Px=350  kgf/m  and 
Kf  =0.3. 

Now  we  will  consider  the  possibility  of  using  the  algorithm 
of  twofold  transfer  for  landers  with  the  designed  ballistic 

parameters:  300  kgf/m2^-  Px  &600  kgf/m2,  0.3$K^.^0.5  in  the  range 

of  entry  velocities  of  5.5  km/s^.VBx^  7.7  km/s. 

During  the  investigation  of  the  accuracy  of  control  for  each 
specific  case  their  own  lines  of  change-over  were  calculated.  The 
results,  based  on  the  appraisal  of  the  accuracy  of  control  depend- 
ing on  load  on  the  frontal  surface  of  the  lander,  lift-drag  ratio, 
and  entry  velocity,  are  given  in  Figure  4.  It  is  evident  that  an 
increase  in  the  load  on  the  frontal  surface  has  practically  no  in- 
fluence on  the  accuracy  of  control,  and  an  increase  in  the  lift-drag 
ratio  worsens  it  somewhat.  Thus,  for  example,  with  an  increase 
of  Px  from  350  to  550  kgf/m2  (K^>  =0.3;  VBX=5.8  km/s)  AVk  changes 
from  45  to  50  m/s,  and  with  an  increase  of  K^.  from  0.3  to  0.5 

8. 
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(Px=350  kgf/m2,  Vbx=5.8  km/s)  AV^  changes  from  45  to  60  m/s. 

This  is  explained  by  the  fact  that  in  all  cases  the  magnitude  of 
relative  error  in  the  lift-drag  ratio  was  selected  the  same  - 10%. 
With  an  increase  of  quality  there  is  an  increase  in  the  absolute 
magnitude  of  error,  which  leads  to  an  increase  in  the  scattering 
of  final  velocity. 

; 

As  it  should  be  expected,  the  magnitude  of  mathematical  ex- 
pectation of  final  velocity  depends  considerably  on  the  reduced 

i 

load  on  the  frontal  surface,  and  an  increase  of  Px  leads  to  an 
increase  of  M(Vfc).  Thus,  if  for  a lander  with  Px=350  kgf/m^ 
(Ktf-=0.3,  V =5.8  km/s)  M(V  ) = 638  m/s,  then  for  a lander  with 
Px=550  kg/m2  the  magnitude  of  mathematical  expectation  of  final 


velocity  increases  to  842  m/s.  With  an  increase  of  lift-drag  ratio 
the  magnitude  of  M(V  ) initially  drops,  and  then,  beginning  with 
a certain  value  of  Kg-  , is  virtually  unchanged.  Thus,  with  an 
increase  of  K^.  from  0.3  to  0.4  (Px=350  kgf/m,  VB><=5.8  km/s)  M(Vh) 
is  reduced  from  638  to  560  m/s,  and  with  a further  increase  of  K£ 
from  0.4  to  0.5  M(V  ) is  reduced  all  told  by  7 m/s.  It  follows 

K 

from  here  that  an  increase  of  the  lift-drag  ratio  above  0.5  for  a 
lander,  which  is  being  controlled  by  the  roll  attitude,  is  not 
advisable . 

It  Is  evident  from  the  data  given  in  Figure  4 that  an  increase 
of  entry  velocity  has  a weak  influence  of  the  accuracy  of  control. 
Thus,  for  example,  with  an  increase  of  VBx  from  5.8  to  7.5  km/s 
AV  is  changed  from  43  to  35  m/s (for  landers  with  Px=350  km/s2  and 
K/j'=0.3).  It  is  also  evident  that  with  an  increase  of  entry  veloc- 


ity  the  magnitude  of  mathematical  expectation  of  final  velocity 
is  reduced  somewhat.  For  example,  with  an  increase  of  entry  veloc- 
ity from  5.3  to  7.5  km/s  M(VU)  is  reduced  from  638  to  600  m/s. 


Figure  4 . 

Key:  (1)  [m/s];  (2)  [kgf/m2] ; (3)  [km/s]. 

Thus  the  materials  presented  show  that  the  control  algorithm 
which  utilized  the  lines  of  change-over  can  be  used  in  the  construc- 
tion of  systems  for  control  of  final  velocity  for  a wide  class  of 
landers  and  for  different  entry  velocities.  In  this  case  the  con- 
trol algorithm  ensures  an  accuracy  of  no  worse  than  50  m/s,  and 
the  increase  in  final  velocity  in  respect  to  optimal  control  com- 
prises 120  m/s  on  the  average. 

Submitted  19  Oct  1971 
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